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Sedimentation-diffusion equilibrium density profiles of suspensions of charge-stabilized colloids are 
calculated theoretically and by Monte Carlo simulation, both for a one-component model of colloidal 
particles interacting through pairwise screened-Coulomb repulsions and for a three-component model 
of colloids, cations, and anions with unscreened-Coulomb interactions. We focus on a state point 
for which experimental measurements are available [CP. Royall et al., J. Phys.: Cond. Matt. 
17, 2315 (2005)]. Despite the apparently different picture that emerges from the one- and three- 
component model (repelling colloids pushing each other to high altitude in the former, versus a 
self-generated electric field that pushes the colloids up in the latter), we find similar colloidal density 
profiles for both models from theory as well as simulation, thereby suggesting that these pictures 
represent different view points of the same phenomenon. The sedimentation profiles obtained from 
an effective one-component model by MC simulations and theory, together with MC simulations 
of the multi-component primitive model are consistent among themselves, but differ quantitatively 
from the results of a theoretical multi-component description at the Poisson-Boltzmann level. We 
find that for small and moderate colloid charge the Poisson-Boltzmann theory gives profiles in 
excellent agreement with the effective one-component theory if a smaller effective charge is used. 
We attribute this discrepancy to the poor treatment of correlations in the Poisson-Boltzmann theory. 

PACS numbers: 82.70.Dd, 64.60.Cn, 64.10.+h, 82.39.Wj 
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INTRODUCTION 



Colloidal particles with a density different from that of 
the dispersive medium sediment because of the gravita- 
tional force. At fixed temperature T, the resulting non- 
homogeneous equilibrium distribution is a consequence of 
the balance between energy and entropy of the different 
chemical species involved. This equilibrium is character- 
ized by measurable density profiles [3, 0]. In the case 
of sufficiently dilute suspensions those profiles obey the 
barometric law p{x) oc exp(— x/L), with L = hsT/rug 
the gravitational length, m the buoyant mass, g the 
gravitational acceleration, ks the Boltzmann constant, 
T the absolute temperature, and p(x) the number den- 
sity of colloids at altitude x. In dense systems, with non- 
negligible colloidal interactions, strong deviations from 
the barometric law have been observed, e.g. for col- 
loidal hard-spheres at packing fractions up to and be- 
yond the freezing point [l], Q ■ More surprisingly (at least 
initially) were the strong deviations from the baromet- 
ric law in rather dilute suspensions of highly charged 
colloids at low salinity The measured density pro- 
files suggested an extreme enhancement, by at least one 
order of magnitude, of the apparent mass of the col- 
loids 0, HI- This system was theoretically analyzed in 
terms of a three-component model of colloids and mono- 



valent cations and anions, for which Poisson-Boltzmann 
(PB) theory in gravity revealed a self-consistent electric 
field that pushes up the colloids to high altitude against 
gravity, thereby reducing the apparent mass as observed 
experimentally Q . However, another explanation was 
given more recently in Ref. [6i], where hydrostatic equi- 
librium in a one-component system of colloids interacting 
through pairwise screened-Coulomb repulsion was con- 
sidered. In the present paper we investigate the relations 
between these two pictures in more detail by considering 
both models (colloid-cation-anion mixture and colloids- 
only system) by theory as well as Monte Carlo simulation. 
For the one-component approach we use a model based 
on effective pairwise screened-Coulomb interactions. The 
profiles are obtained from the solution of the hydrostatic 
equilibrium equation, that uses the isothermal compress- 
ibility obtained from the solution of the Ornstein-Zernike 
(OZ) equation with the rescaled mean spherical approx- 
imation (RMSA) closure [7]. This approach is similar 
to that introduced in [f|. In this case, entropic and 
electrostatic effects are implicitly included in the struc- 
ture of the suspension (see next section). Furthermore, 
we perform MC simulations for this model. Within the 
multi-component picture, we approach the problem us- 
ing the Poisson-Boltzmann theory introduced in Ref. [f| , 
which explains the non-barometric profiles in terms of a 
macroscopic electric field that appears as a consequence 
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of a charge inhomogeneity. We also performed simula- 
tions of this system using the primitive model in grav- 
ity. These particular simulations require a substantial 
amount of CPU time since a considerable number of 
micro-ions has to be taken into account to mimic the sub- 
stantial salt concentration in the suspension. We focus on 
a state point for which experimental information is avail- 
able, and compare the profiles obtained from the men- 
tioned theoretical and numerical approaches with pub- 
lished measurements Q. 

The paper is organized as follows. In section 2 we in- 
troduce the one-component model, the structure factor of 
the suspension determined through RMSA, and we dis- 
cuss some details regarding the simulation technique of 
the one-component model. In section 3 we briefly revise 
the PB theory for sedimentation, introduce the primitive 
model in gravity and some aspects regarding the simula- 
tions. Results and discussion are presented in section 4 
and 5 respectively, where we compare the different theo- 
retical and numerical results with experimental data. A 
summary and conclusions are gathered in section 6. 



EFFECTIVE ONE-COMPONENT MODEL 

Let us consider a system consisting of charged spheres 
of diameter a = 2a, mass M, and electric charge — Ze, 
in osmotic contact with a reservoir of 1:1 electrolyte 
with salt concentration 2p s . The solvent has mass den- 
sity pi and is characterized by an electric permittivity 
e. Let us assume also that the dielectric constant of the 
spheres and the electrolyte are identical to avoid electro- 
static image effects and Van der Waals forces between 
the spheres. Assuming pairwise effective colloidal inter- 
actions, the Hamiltonian of the effective one-component 
system of colloids in the presence of gravity is given by 



H- 



N 



N 



mgx. 



(1) 



i<j 



where the first term in the right-hand side is the potential 
energy of colloid i at height Xi measured from the bottom 
of the sample. Here m = M — piira 3 /6 the buoyant 
mass of the colloidal particles, and v(r) is the familiar 
screened-Coulomb potential 

{oo , if r < a, 

Z 2 exp( K q) X B (2) 
— — r~ exp(-KT), if r>a, 
(1 + nay r 

with (3 — (ksT)^ 1 where ks is the Boltzmann con- 
stant and where Xb — ^f- is the Bjerrum length, k = 
\/8tt\bPs is the inverse screening length, and r is the 
distance between centers of colloidal particles. Under 
isothermal conditions and for small density gradients, the 
osmotic pressure of the suspension depends only on the 



local number density of colloids p(x). The latter is de- 
termined from the non-linear differential equation that 
follows from inserting p(x) into the hydrostatic equilib- 
rium equation dli{x)/dx — —mgp(x) with II the osmotic 
pressure of the suspension with respect to the salt reser- 
voir. This yields 



dp{x) Xt{p{x)) 
dx L 



p(x) 
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where \ T X 



\~§p) * s ^ ne isothermal compressibility 
of the bulk fluid and L is the gravitational length defined 
above. The sedimentation profiles can be obtained by 
solving ([3]) if the function xt(p) is known for the relevant 
density regime. In order to determine xt{p) we use the 
well-known Kirkwood-Buff relation xt = lim S(q) with 

S(q) the structure factor as calculated by Hansen, Hayter 
and Penfold within the RMSA closure of the Ornstein- 
Zernike equation 0, [2l[ ■ By this procedure the sedimen- 
tation profiles are determined solely from the structure of 
the effective one-component bulk fluid. Notice that such 
a scheme was applied successfully to explain the mea- 
sured hard-sphere density profiles in Ref. [l[. For later 
comparison we also consider an alternative expression for 
Xt{p)i that is based on the Donnan equation of state as 
e.g. given in Ref. [5(. This yields 



Xt 
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which features the high-density or low-salt limit \t ip) — 
1 + Z for Zp 3> 2p s , such that insertion into Q yields 
p(x) tx exp[— x/(Z + 1)L], i.e. an effective gravitational 
length that is a factor Z + l larger than that in the baro- 
metric law [5[. The remaining task in order to find the 
sedimentation profiles is to insert xt into ([3]) and to solve 
the non-linear equation numerically on an x-grid. 

In addition we perform standard Monte Carlo simula- 
tions of a system described by the interaction Hamilto- 
nian |T]) for the parameters Z = 76, colloid diameter 
a = 1.91^tm, Bjerrum length Xb = 10.4 nm, screen- 
ing parameter na — 1.2 and average colloidal packing 
fraction fj = H^ 1 j£* rj(x)dx — 0.0053 with the height 
H = 50.92cr. The experimental screening parameter sat- 
isfies kg — 1.2. These parameters are identical to those of 
the experimental system studied in Ref. Q , where Z = 76 
stems from the best fit of the experimental density pro- 
file with a theoretical prediction based on the primitive 
model (see below). The dimensions of the rectangular 
simulation box are lOcr x 10a x H. We checked that the 
horizontal area was large enough to exclude finite-size ef- 
fects. We employed periodic boundary conditions in the 
horizontal directions, in the vertical directions the sys- 
tem is bounded by hard walls that exclude the centers of 
colloids at x < and x > H. 
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THE PRIMITIVE MODEL IN GRAVITY 

As mentioned in the introduction, a different approach 
to study sedimentation profiles is to consider each chem- 
ical species separately, namely colloids (c), coions (-) and 
counterions (+). The Hamiltonian of the system can be 
written as 



H= H r 



Ha + H c 



N 

^mgx t 

i=l 



(5) 



where the first three terms in the right-hand side include 
colloid-colloid, ion-ion and colloid-ion pairwise interac- 
tions, respectively, and the last term is the gravitational 
energy of the colloids introduced in equation fl} ; the ions 
are assumed to be massless. In this model, the electro- 
static pair interactions are of the form 



oo , 



if r < ov 



ZiZj^f, if r>a tj , 



(6) 



with CTfc and Zk the diameter and the charge number of 
species k = {c, +,— }, i.e. Z c = — Z,Z + = 1,Z- = —1, 
a c — a and <7 + = er_ <C a. The number of particles 
are denoted by N k , i.e N c = N and N + = N- + ZN 
for charge neutrality reasons. This three-component 
model can be studied within Poisson-Boltzmann theory 
which relates the density profiles p(x),p±(x),p-(x) of 
the colloids cations and anions respectively, to the local 
electrostatic Donnan potential ip(x) through 

P±(x) = P s exp[±<t>{x)]; (7) 
p(x) = po exp[-x/ L + Zcj>(x)] ; (8) 
tf>"(x) = K 2 smh4>(x) +4tt\ b Zp(x), (9) 

with 4>(x) defined by the dimensionless combination 
4>{x) = eip(x) l 'ksT and where po is a normalization con- 
stant. Here a prime denotes a derivative w.r.t. x. Un- 
der appropriate conditions, typically Z 2 p(x) ^> 2p s , it 
was found that <j){x) is a linear function of x in macro- 
scopically large parts of the system, i.e. there is a con- 
stant electric field that lifts the colloids to higher al- 
titudes than expected on the basis of their mass 
This result stems both from numerical solutions of 0- 
© and from the " Laplace-Boltzmann" equation, where 
© is replaced by the local charge neutrality condition 
k 2 sinh0(x) + AttXbZp(x) — 0. Note that by combining 
the latter equation with (0 and (|8|) one recovers the hy- 
drostatic equilibrium condition (|3|) with \t given by (|4|). 
On the other hand, the set of equations I©-© can be 
solved numerically in order to determine the local electro- 
static potential 4>{x) together with the equilibrium pro- 
file p(x) by an iterative procedure as pointed out in 
Such a procedure requires two boundary conditions, e.g. 
(f>'(0) = 4>'{H) = 0, where H is the height of the solvent 
meniscus. 



A system described by the Hamiltonian ([5]) was simu- 
lated in a rectangular box of horizontal area 9a x 9a and 
height H — 50<t. The vertical coordinate x is restricted 
to x £ [0, H], and periodic boundary conditions are only 
applied in the horizontal plane and not in the vertical 
direction. In order to be as closely as possible to the ex- 
periments of Ref.0], we considered colloids with charge 
Z = 76, diameter a = 1.910/im, gravitational length 
L = 2.41(7, Bjerrum length = 10.4nm, and average 
colloidal packing fraction rj = H^ 1 J^ 1 rj(x)dx = 0.0053 
(with H = 50a). The experimental screening parame- 
ter satisfies na = 1.2. This state point is realized, with 
the present box size and shape, by the number of col- 
loids -/V c =12 and the number of positive and negative 
ions N + =13516 and iV_ =12604, respectively. In order 
to take the long-range electrostatic interactions into ac- 
count we have employed a combination of Ewald summa- 
tion in a slab geometry with the lattice method proposed 
by Panagiotopoulos and Kumar The parameters of 
the Ewald summation and lattice method are the same 
as those in reference [13] ■ Note that the large number 
of ions, which have to be included to represent the low 
but yet substantial screening parameter, makes the sim- 
ulations extremely time consuming. A typical simula- 
tion consists of 10 5 MC cycles. A cycle consists of 0.9Af 
trials to move a randomly chosen colloid and 0.1 N tri- 
als to move a randomly chosen colloid or microion, with 
J\f = N + N + + N- the total number of microions and 
colloids in the system. In a dense system of microions, 
a simple Monte Carlo move of a colloid would almost 
certainly result in an hard-core overlap with one of the 
microions. In order to avoid such overlaps we use a clus- 
ter move technique, where ions that overlap with the new 
colloid position are moved into the space left empty by 
the displaced colloid, more details on this technique can 
be found in Ref. 18, [2(|. The percentage of accepted 



moves of each component (colloids and microions) was 
maintained at about 40%. To check if the system was 
equilibrated, the average altitude of the centers of mass 
of the colloids was monitored in the simulation; when 
the center of mass was not stable, further equilibration 
was performed before taking measurements. A final sim- 
ulation with 2 • 10 5 MC cycles was performed to obtain 
averages. 



RESULTS 

Figure 1 shows a first comparison of theoretical predic- 
tions based on the one-component and three-component 
model. We see several theoretical sedimentation profiles 
as a function of the altitude x, corresponding to colloids 
of diameter a = 1.91/xm, and three different salt concen- 
trations characterized by screening parameters na = 0.8, 
1.2, and 1.6. All profiles shown in Fig.l are for the same 
gravitational length L = 2.41<r, average packing frac- 
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FIG. 1: Colloidal sedimentation profiles based on hydrostatic 
equilibrium © calculated using the RMSA-based compress- 
ibility of the one-component Yukawa model, compared to 
those based on the multi-component PB theory ©-©, for 
the screening parameters no — 0.8, 1.2 and 1.6. The col- 
loidal charge is Z — 76, the colloid diameter a = 1.91^im, 
the Bjerrum length is As = 10.4 nm, the gravitational length 
L — 2.41<t, the average colloid packing fraction is rj — 0.0053, 
and the sample height is H = 50<r, as reported in Ref. 
The dashed curve is the barometric distribution with the same 
normalization. The inset shows the RMSA-based compress- 
ibility. 



tion rj — 0.0053, sample height H — 50a, and Bjerrum 
length Xb = 10.4nm. For each k<j, the colloid density 
profile is calculated for both the effective one-component 
model based on the solution of the hydrostatic equilib- 
rium equation ([3]) using the isothermal compressibility 
obtained from the RMSA closure, as well as from the 
multi-component PB theory described by equations |7J)- 
([9]). We also show, for the sake of comparison, the cor- 
responding barometric profile obtained from |8]) and © 
in the case of uncharged colloids (Z = 0), for the same 
normalization. The inset shows the compressibilities as 
a function of the colloid density, as obtained from the 
solution of the OZ equation within the RMSA closure. 
At zero density all the compressibility curves reduce to 
the ideal-gas compressibility, and with increasing col- 
loid density the electrostatic repulsions manifest them- 
selves as a reduction of xt- the weaker the screening, 
the stronger the effective colloidal interactions. We note 
that each one-component Yukawa system yields steeper 
density profiles than those of the corresponding three- 
component model, for all na considered here, i.e. the 
one-component systems have a relatively small average 
altitude and a relatively low density at higher altitudes. 
We will argue in more detail below that the source of the 
difference between the one- and three-component predic- 
tions is mainly due to the poor representation of the 
colloid-colloid correlations in the three-component PB 
theory. 
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FIG. 2: Colloidal sedimentation density profiles for colloidal 
charge Z — 76 stemming from multi-component PB theory 
0-©; from the Donnan compressibility © combined with 
hydrostatic equilibrium ©, and simulations of the Primitive 
Model in gravity, for the parameters of Fig.l and no — 1.2. 
The symbols denote the experimental measurements from 
Ref.0]. 



Figure 2 shows experimentally measured density pro- 
file of Ref. Q compared to density profiles as obtained 
from the multi-component models: the three-component 
PB theory of equations CO)-© and simulation of the 
primitive model in gravity as introduced in section 3. 
In spite of the fact that the primitive model simulation 
was equilibrated during about one year CPU time and 
the measurements were performed over four months CPU 
time, the level of noise in the raw data is still quite high. 
Therefore also a smoothed curve of the simulated pro- 
file is shown to facilitate comparisons. The difficulty in 
obtaining good statistics in this particular simulation is 
due to the fact that the total number of colloids in the 
system is exceedingly small (N c = 12), whereas the total 
number of particles in the system is rather large (26132), 
most of them salt ions needed to achieve the required 
screening parameter condition. That the Donnan-based 
density profile is accurate when compared to the experi- 
ments is only due to the fact that the experimental value 
Z = 76 [33| stems from a fit to PB theory Q, which is 
equivalent to the Donnan equation of state in the local 
neutrality " Laplace-Boltzmann" limit as explained below 
Eqs.®-©. In other words, Z = 76 is close to a best fit 
to the Donnan equation of state. 

In figure 3 we see a first comparison of the sedimen- 
tation profile obtained experimentally with sedimenta- 
tion profiles calculated using the effective one-component 
models: simulation of the Yukawa system and the RMSA 
approach of section 2. We also include the profiles ob- 
tained from the multi-component PB theory for the sake 
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FIG. 3: Sedimentation profiles of the effective one-component 
Yukawa model calculated using both RMSA and standard MC 
simulations, compared with the experimental measurements 
of Ref. Q] and the three-component PB theory, for Z = 76 
and all other parameters as in Fig. 2. Note the close agree- 
ment between the two Yukawa results, and their small but 
systematic deviation from the experiments and the PB the- 
ory. 

of comparison. The contrast between the simulations and 
the experimental curve seems to reveal a systematic de- 
viation such that the simulated and RMSA profiles are 
actually somewhat too steep. Indeed, when we allow Z 
to be a fit parameter in the one-component Yukawa sys- 
tem, keeping all the other parameters equal, it turns out 
that Z = 94 gives best agreement of the one-component 
models with the experimental profile. It is tempting to 
conclude, therefore, that Z = 76 gives merely a best fit 
to the experiment within PB theory given by Eqs.((7|)- 
©, which (to a large extent) ignores colloid-colloid cor- 
relations, whereas inclusion of these correlations (as in 
the simulations of the primitive model and that of the 
Yukawa system, and in the RMSA-based theory) gives 
rise to a density profile that is systematically steeper in 
comparison with the experiment. 

In figure 4 we show the resulting sedimentation pro- 
files based on the Yukawa potential simulations and the 
RMSA closure for Z = 94. For comparison we also plot 
the multi component PB model for Z = 94 revealing a 
relatively poor agreement with the other curves. In this 
case the PB approach clearly fails to reproduce quantita- 
tively the sedimentation profiles. On the other hand, the 
agreement of the experimental profile with the effective 
one-component models is good, except that as mentioned 
before, the simulated profile exhibits much more struc- 
ture close to the hard-wall that represents the bottom 
of the sample in the simulations — this packing effect is 
not captured by the local-density approximation that un- 
derlies the hydrostatic equilibrium condition, and is not 
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FIG. 4: As Fig.3, but with Z = 94, such that the Yukawa 
model profiles (both RMSA and simulated) fit the experi- 
mental profile best. The PB profile is clearly less accurate 
now. 

seen in the experiment because the actual sample extends 
beyond the plotted x-range. 

In Fig. 5 we show sedimentation profiles as obtained by 
simulations of the primitive model and of the Yukawa 
model, compared with those by the RMSA approach, all 
for Z = 76. The agreement is perhaps a bit less quanti- 
tative than one would have expected. One of the reasons 
that the density in the primitive model is considerably 
higher in i/cr e (10,25) is due to the structure close to 
the hard wall at the bottom near x — 0, as shown in the 
inset of Fig. 5, where the two Yukawa systems reveal a 
larger net adsorption than that of the primitive model, 
albeit for different reasons: the simulated Yukawa system 
shows a strong peak at x = a/2 while the RMSA-based 
profile continues to be nonzero down to x = 0. Given 
that we imposed that fj is identical in all cases, there 
must also be a region in space where the density in the 
primitive model exceeds the other two; the order of mag- 
nitude of the integrated differences over x/cr G (10, 25) is 
indeed similar to the negative of that over x/a £ (0,3). 
Another reason for these differences might be the poor 
statistics and slow equilibration of the primitive model 
simulations. Recall that the present data are based on 
about one year of CPU time, so considerable extensions 
and more checks are not easily obtained. This also pre- 
vented us from performing primitive model simulations 
for Z = 94. 



DISCUSSION 

The first observation from Fig. 3 should be the gross 
agreement between the experiments and all calculated 
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FIG. 5: Sedimentation density profiles as obtained from sim- 
ulations of the Primitive Model in gravity and simulations of 
the Yukawa fluid compared with the one-component RMSA 
model for colloidal charge Z = 76 and all other parameters 
as in Fig. 2. The difference between the profiles can partly be 
attributed to the structural differences close to the bottom 
at x — as shown in the inset, where the simulations re- 
veal hard- wall induced structure that is not captured by the 
RMSA-based theory, and perhaps partly by slow equilibra- 
tion and poor statistics in the simulations due to the small 
number of colloids. 



and simulated profiles. In all cases we have na = 1.2, 
Bjerrum length = 10.4 nm, colloid diameter a = 
1.91/xm, and gravitational length L = 2 Ala. The col- 
loidal charge is taken as Z — 76, and the measured 
packing fractions are in the range < r/(x) < 0.02, 
where rj(x) — (ir/6)a 3 p(x). This regime is such that 
Zp(x)/2p s < 0.17 for all x, i.e. even at the highest den- 
sity the ion concentration is dominated by the reservoir 
salt concentration 2p s , such that the screening constant is 
indeed essentially a constant independent of the height or 
density, as implicitly assumed in equation @. A closer 
look, however, shows that even though the simulations 
and the RMSA result of the Yukawa system are very close 
to each other (except at the bottom where packing effects 
affect the simulations), they both deviate systematically 
from the experiment: the former two are too steep and 
have too low a density at higher altitudes. From the fact 
that the RMSA-based profile and that of the Yukawa 
simulations are so close to each other, one could conclude 
that they are mutually consistent and both accurate, and 
that their deviation from the experiment is mainly due 
to the present choice of Z = 76, which was based on the 
fitting to the PB theory of Eqs.tfT])-©. This fitting is not 
optimal due to the inadequacy of the present PB theory 
to account for correlations among the different species in 
the system. In particular, PB theory overestimates the 
colloid density at high altitudes for a given value of Z. 



This happens for charges Z > 50, whereas for smaller 
values of the colloid charge the agreement between the 
profiles obtained from the PB and RMSA approaches is 
excellent, as far as an effective (smaller) value for the 
charge is used in the PB approach as discussed in detail 
below. For the parameters of present interest, fitting the 
experimental density profile to that of the Yukawa sys- 
tem treated within the RMSA closure, we concluded that 
Z = 94 gives the best fit. 

Given that a Yukawa system with a colloidal charge 
Z = Zrmsa = 94 is best described within PB theory by 
a colloidal charge Z = Zpb = 76, for the present system 
parameters, it is interesting to investigate the relation 
between Zrmsa and Zpb for other values of the colloid 
charge. In figure 6 we plot the ratio Zpb /Zrmsa for 
< Zrmsa < 200, for three screening constants while 
all the other parameters are left unchanged. Figure 7 
shows the corresponding density profiles for no = 1.2 and 
Zrmsa — 200, 94, and 40. One can conclude that for 
Zrmsa ^ 50, PB theory reproduces the RMSA- Yukawa 
sedimentation profiles very accurately provided the col- 
loidal charge Zpb is reduced by up to 20% of Zrmsa- 
Note that the average density here is low enough that 
the profile in the limit of Zrmsa — * becomes essen- 
tially barometric; at higher average packing fractions one 
also expects Zpb I Zrmsa 7^ 1 in this limit, due to hard- 
core effects that are not accounted for properly in the 
PB theory. The required reduction of Zpb / Zrmsa with 
increasing Zrmsa > 50 exceeds 20%, and, in addition, 
the quality of the best-fitting PB profile (quantified by 
the mean-squared deviation) becomes slightly less satis- 
factory as is reflected by the increase of the error bars 
in figure 6 with increasing Zrmsa, this is also shown in 
figure 7. 

The difference between Zrmsa and Zpb in the present 
system is of course considerable and significant, but not 
qualitative. The seemingly different mechanisms that un- 
derly the lifting of the colloids to higher altitudes than 
given by the barometric distribution, as predicted by the 
three- and one-component theory, should therefore be ac- 
tually equivalent: the self-consistent electric field that 
is generated by a net charge-imbalance at the bound- 
aries of the 3-component system such that the colloids 
are pushed upwards is merely another way of describ- 
ing a pairwise screened- Coulomb repulsion that pushes 
the colloids apart to higher altitudes in a one-component 
model. This is in line with conclusions in Ref.[6j. 

We have attributed the difference between the best 
fit for Z based on Poisson-Boltzmann theory and the 
other three methods (simulations of the primitive model, 
and simulations and RMSA theory of the one-component 
Yukawa system) to the poor account of correlations in the 
Poisson-Boltzmann theory. In principle, however, there 
could be other sources that cause such a difference, e.g. 
charge renormalization and hard-core exclusion effects for 
the screening ions. Charge renormalization due to nonlin- 
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FIG. 6: Ratio of the best-fitting colloidal charge Z = Zpb 
(see text) and that of the RMSA charge Zrmsa, as a function 
Zrmsa G (0,200), with a, As, L, fj, and H as in Fig.l, for 
various screening parameters kg. Note that PB theory is 
increasingly better for lower colloidal charges and lower salt 
concentrations. The lines are mere guides to the eye. 
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FIG. 7: Sedimentation profiles obtained from RMSA theory 
compared to those of PB theory using the best-fitting Zpb 
from Fig. 6, for kg — 1.2 and all the other parameters as in 
Fig.6. Note that the quality of the fit detoriates slightly with 
increasing Zrmsa- 



ear screening effects |26[ is, however, not a candidate here 
to explain the difference for at least two reasons: (i) The 
actual bare charge is usually larger than the renormal- 
ized charge that appears in the prefactor of the screened- 
Coulomb interactions, whereas here the former seems to 
be smaller, (ii) The present parameters here are such that 
ZXg/a < 1, whereas charge renormalization is only sub- 
stantial if this dimensionless combination exceeds about 5 
or so [25, 2(| 29[. A mechanism whereby the effective col- 
loidal charge is increased was discussed in Refs. 



27 



and is based on the hard-core exclusion of the screen- 
ing ions at sufficiently high colloid packing fractions: the 
screening is therefore less effective, which appears as an 
increase of the effective colloidal charge. However, apply- 
ing the analysis of Ref. [28| to the present case gives only 
marginally larger values of the effective charge, by less 
than 1% at the highest density r\ ~ 0.02. In other words, 
it appears that this effect cannot explain the difference 
between Z = 76 and Z = 94, leaving the poor account 
of colloid-colloid correlations in the Poisson-Boltzmann 
theory as the most plausible source of the difference. 

It is interesting to inquire whether the one- and three- 
component models would also produce essentially the 
same sedimentation profiles for other sets of parameters 
than considered here, and whether the hydrostatic equi- 
librium condition J3J and the Poisson-Boltzmann theory 
([I])-© for the one-component and the three-component 
case, respectively, produce reliable profiles in all circum- 
stances. In order to answer these questions we consider 



the regimes of extremely low and extremely high salt con- 
centrations. First, the hamiltonian JT]) for the effective 
one-component system is pairwise additive, which is ex- 
pected to be a good approximation for the present pa- 
rameter set where no~ > 1 and Zpj1p s < 1, i.e. the 
range of the interactions is smaller than the size of the 
particles and the background electrolyte concentration 
dominates the counterion concentration. At lower salt 
concentrations, such that w <C 1 and/or Zp/2p s > 1, 
one would expect effective many-body interactions to be- 
come relevant 14, 2j|, such that ([1]) is not necessarily a 
reliable effective hamiltonian anymore. In such an ex- 
tremely low-salt regime, which is realized in salt free sys- 
tems, the Poisson-Boltzmann theory proved to be quanti- 
tatively accurate, at least in comparison with simulations 
liL 32 ] at low Coulomb couplings. It is interesting to see 
if the pairwise one-component description is capable of 
describing the density profiles in this case. We wish to 
stress here that the possible break-down of the pairwise 
screened-Coulomb picture does not imply that the sys- 
tem can no longer be seen as a one-component system in 
hydrostatic equilibrium as described by ((3j with a com- 
pressibility that follows from the Kirkwood-Buff relation 
S(0) = xt- these relations remain valid (the former only 
within the local density approximation, but given the 
long screening length in the extremely low-salt regime 
this approximation is probably accurate). The break- 
down would "merely" imply that it is not obvious how 
to calculate the compressibility or the structure factor 
without detailed knowledge of the effective hamiltonian. 
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Second, let us consider the opposite high-salt regime such 
that no- 3> 1 and Zp <C 2p s . In this regime the elec- 
trostatic interactions are completely screened over dis- 
tances much smaller than the colloidal diameter, such 
that the effective one-component system is essentially a 
(pairwise) hard-sphere system (for water at room tem- 
perature at least, where ion-ion correlations are not all 
that important). In this regime the one-component de- 
scription based on (|3|) is far superior over the PB theory 
of Eqs. ([7|)-(l9]). This is directly seen by regarding the 
Z = limit of equations ([T])-©, which reduce to 4>(x) = 
and p{x) = pa exp[— x/L], i.e. the sedimentation pro- 
files become barometric; the hard-core correlations are 
left-out completely from this theory. By contrast, the 
RMSA closure is, in this hard-core limit, equivalent to the 
Percus-Yevick closure, and in combination with © the 
density profiles of hard-sphere sedimentation equilibrium 
are well-described P, Moreover, also in the present 
regime with na ~ 1 the one-component theory performs 
better. We are currently working on the formulation of 
a theory that is able to describe sedimentation density 
profiles in both the high-salt and the low-salt regime on 
the same footing. 

SUMMARY AND CONCLUSIONS 



seen that PB theory actually overestimates the colloid 
density at high altitudes compared to the corresponding 
Yukawa system, for identical values of Z , at the parame- 
ters of interest here. Agreement between PB theory and 
Yukawa systems can be obtained by reducing the col- 
loidal charge in the PB theory compared to that of the 
Yukawa model. For small values of the colloid charge, 
Z < 50 or so, the agreement between the resulting pro- 
files obtained from the PB and RMSA approach is truly 
excellent, for larger charges up to say Z = 200 the agree- 
ment is still good though somewhat less quantitative as 
regards the functional form. The CPU time required for 
the simulation of the multi-component primitive model 
and the effective one-component Yukawa model varies 
between about one year in the former case and one hour 
in the latter. This shows that theory and simulations of 
sedimentation profiles on the basis of the effective one- 
component potential and the Poisson-Boltzmann theory 
(possibly with a reduced effective charge when colloid- 
colloid correlations are important) are considerably more 
efficient than primitive model simulations. In spite of the 
fact that we have only considered a particular case study, 
this is presumably true in general, with the Yukawa 
model probably more accurate when na > 1 while PB 
theory could be more accurate or efficient when na <C 1. 
This will be investigated in more detail in a future study. 



In this paper we have studied sedimentation equilib- 
rium of charge-stabilized colloids at non-zero salt con- 
centration. We compared experimental results with the- 
oretical and simulated profiles obtained on the basis of 
two models. On one hand a multi-component model of 
point-like colloids, cations, and anions interacting with 
bare Coulomb potentials. For this model we considered 
a Poisson-Boltzmann theory of the three-component mix- 
ture and performed MC simulations using 12 colloids 
and a total of about 26132 particles to guarantee the 
electroneutrality of the system. On the other hand we 
considered an effective one-component model of colloids 
interacting by an effective screened-Coulomb potential. 
For this model we employed a theory based on hydro- 
static equilibrium, where the isothermal compressibility 
is given by the Kirkwood-Buff relation as obtained from 
the solution of the Ornstein-Zernike equation with the 
rescaled mean spherical approximation (RMSA) closure 
for the screened-Coulomb potential. For the effective 
one-component Yukawa model, we also performed sim- 
ulations of sedimentation profiles. 

The sedimentation profiles obtained from the one- 
component RMSA theory, simulations of the Yukawa sys- 
tem, and simulations of the primitive model are essen- 
tially consistent among themselves but differ from the 
results of the Poisson Boltzmann theory. The PB theory 
shows good agreement with the experiments only because 
the numerical value of the charge was estimated as to give 
the best fitting according to this theory. In fact, we have 
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